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Evolution of inflation-generated magnetic field through phase transitions 
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We study the evolution of an inflation-generated magnetic field, due to its coupling to fluid mo- 
tions, during cosmological phase transitions. We find that the magnetic field stays almost unchanged 
on large scales, while on small scales the spectrum is modified in such a way that power at small 
scales becomes progressively suppressed. We also show that the magnetic field generates turbulent 
motions in the initially turbulence-free plasma. On large scales, the slope of the resulting kinetic 
energy spectrum is consistent with that of white noise. 

PACS numbers: 98.70.Vc, 98.80.-k 



I. INTRODUCTION 



The origin of the coherent large-scale (~ 10 kpc) part of 
galactic magnetic fields, of /xG strength, is under active 
discussion [ll-Q. On larger, Mpc scales, until recently 
there were only upper limits, the most restrictive being of 
order a few nG depending on the observational technique 
used to measure the intergalactic magnetic field strength 
Recently, there have been a number of published 
lower limits on a putative large-scale magnetic field of 
strength HT 1 * 1 fG (1 fG = 1CT 15 G) 0, or possibly two 
orders of magnitude smaller [7J. 1 

Almost certainly, the galactic fields are the amplified 
remnants of significantly weaker "seed" magnetic fields. 
Quantum mechanical fluctuations during inflation [9( is 
a leading candidate for generating the needed seed mag- 
netic field To generate a large enough seed mag- 
netic field through quantum mechanical fluctuations dur- 
ing inflation, conformal invariance must be broken dur- 
ing inflation. A simple, realistic, illustrative model cou- 
ples the abelian vector field with field strength tensor 
F^ v to the scalar inflaton field <f> through a dilaton-like 
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1 These techniques for limiting a large scale cosmological magnetic 
field might be unreliable [§[, but see their Sec. 4 where they note 
that more work will be needed to firm up these arguments and 
to determine whether the techniques used to establish the lower 
limits arc indeed unreliable. 



coupling, generalizing the Maxwell lagrangian density to 
e^F^F^ where a is a parameter [111, ll2( . In the case 
of power-law inflation, and depending on the value of a, 
this can result in a large enough seed magnetic field to 
explain the observed galactic magnetic fields. This is an 
observationally viable model. For a more detailed de- 
scription of the model see Sec. II below. 

After the end of inflation, such an inflation-generated 
magnetic field will be correlated over super-Hubble- 
radius scales. It would induce observable signatures 
in the cosmic microwave background (CMB) radiation 
anisotropies at the epoch of recombination (the last scat- 
tering surface) if its current amplitude on Mpc scales is 
of the order of a nG [HJ]. 2 The properties of an inflation- 
generated primordial seed magnetic field depend on the 
parameters of the inflation model. If cosmological ob- 
servations confirm the presence of an inflation-generated 
magnetic field, these measurements could be used to 
probe the physical conditions during inflation, including 
the shape of the inflaton potential energy density as well 
as the coupling between the inflaton and the vector gauge 
field. 

To check the consistency of the model, the primor- 
dial magnetic field shape and amplitude should be mea- 
sured in as many ways as possible. The simplest limit 
arises from the cosmological expansion dynamics during 
big bang nucleosynthesis. This requires that the energy 
density of the magnetic field should not be larger than 



2 The effects of an homogeneous magnetic field on the CMB 
anisotropy, and the resulting non-Gaussianity, are discussed in 
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about 10% of the radiation energy density. This limits 
the present (inflation-generated) magnetic field strength 
to less than a few /iG, if the primordial magnetic field 
was generated prior to or during big bang nucleosynthe- 
sis, and was not damped or amplified by a magnetohy- 
drodynamic (MHD) or some other process and so stays 
frozen into the plasma [l5j . 

In addition to the CMB temperature anisotropies 
that a primordial magnetic field induces (as mentioned 
above), such a field will Faraday-rotate the CMB po- 
larization anisotropies (16|. Currently available Faraday 
rotation data give a bound on the primordial magnetic 
field strength of less than a few nG (for a scale-invariant 
or homogeneous primordial magnetic field). 

Another interesting signature of a cosmological mag- 
netic field is the relic gravitational wave signal generated 
by the anisotropic magnetic stress [13]. The amplitude 
of the induced gravitational waves is determined by the 
magnetic field energy density, so a direct measurement 
of the resulting gravitational wave signal can lead to an 
independent limit on the magnetic field strength; see Ref. 
[l8l | and references therein. 

After the Universe reheats at the end of inflation, the 
plasma that was created then has large conductivity and 
it is conventional to assume that this remains the case 
as the Universe evolves to the present. In this case 
the large-scale cosmological magnetic field behaves as a 
frozcn-in field with an evolution determined by the sim- 
ple, flux-conservation, dilution of magnetic field lines, 
B(x, t) oc B (x)/a 2 (i), where t is the physical cosmic 
time and a{t) is the cosmological scale factor. On the 
other hand, the evolution of a primordial magnetic field 
is a complex process influenced by MHD as well as by 
the dynamics of the Universe [19-22]. In particular, the 
presence of a magnetic field can dramatically affect pri- 
mordial turbulence (e.g., when the turbulence is associ- 
ated with cosmological phase transition bubble motions) 
fl9l I23L H3]. Furthermore, the presence of a magnetic 
field itself might lead to the development of turbulent 
motions, and so affect the turbulence [25l[26j. 

In a recent examination (2?J of the effects of the MHD 
coupling between a primordial magnetic field and turbu- 
lence during a cosmological phase transition, we consid- 
ered two different initial shapes for the spectrum of the 
primordial magnetic field: a single-scale magnetic field 
and a magnetic field with a Batchelor spectrum at large 
scales. In this paper we present a similar analysis for 
modified initial conditions for an inflation-generated pri- 
mordial magnetic field [ll|, coupled via the usual MHD 
equations with the fluid, during the electroweak or QCD 
phase transitions. We consider both non-helical and he- 
lical magnetic field cases. 

We assume that the phase transition bubbles induce a 
typical length scale at which the magnetic field starts 
to interact with the phase transition fluid. The rele- 
vant difference between the electroweak and QCD phase 
transitions is encoded in the difference between values 
of parameters such as the temperature the number 



of relativistic degrees of freedom g+, the bubble number, 
and bubble sizes. We assume initial absence of primor- 
dial turbulence, i.e., we assume that the plasma is ini- 
tially at rest (although it is possible to generate turbulent 
motions through bubble collisions and nucleation 28]). 
The characteristic parameter of the primordial magnetic 
field is the r.m.s. Alfven velocity va — B / ^/167rp rac j/3. 
Here, p ra( j — /^thermal is the radiation energy density. 
We use Slrad^o = 2-56 x 10~ 5 , where Sl ra d is the ra- 
diation energy density parameter and h is the Hub- 
ble constant in units of 100 km s" 1 Mpc -1 , for a cur- 
rent CMB temperature Tq = 2.74 K. At temperature 
T*, PradCT*) = 7r 2 3*(T*) 4 /30. The Alfven velocity does 

not depend on T* but is weakly dependent on g+, i.e., 

-1/6 
va oc 3* 

In our previous simulations [27j we studied phase- 
transition-generated magnetic fields coupled to a rela- 
tivistic fluid and discovered that equipartition between 
kinetic and magnetic energy densities is reached within 
reasonably short times. In this paper our main purpose 
is to consider different initial conditions. In particular, 
in the case when the magnetic field is generated during 
inflation, we investigate the kinds of turbulent motions 
that result from the coupling of the magnetic field with 
the fluid, and determine how this affects the evolution of 
the field itself. We show that the presence of a magnetic 
field on large scales ensures a rapid rise of the velocity 
field on large scales. On the other hand, magnetic field 
decay on large scales occurs at slow rates. 

The structure of our paper is as follows. In Sec. II 
we briefly describe magnetic field generation during in- 
flation. In Sec. Ill we discuss the phenomenological cou- 
pling of the magnetic field to the turbulent plasma. In 
Sec. IV we present our numerical simulation results. We 
discuss our results and conclude in Sec. V. 



II. INFLATION-GENERATED MAGNETIC 
FIELD 

A. Non-helical magnetic field 

Any field included in the lagrangian density during in- 
flation will be produced by quantum-mechanical fluctu- 
ations. These fluctuations are then stretched by the cos- 
mological expansion, leaving the Hubble radius during 
inflation and re-entering it at much later times. 

Adding the standard abelian vector field lagrangian 
density F^F^ / e 2 (where e is the vector potential cou- 
pling constant) to the usual general relativity and scalar 
inflaton field lagrangian density results in electric and 
magnetic fields being generated during inflation. Of 
course, these are not the usual low-energy electric and 
magnetic fields; rather, at the very least, they are hyper- 
charge electric and magnetic fields. During inflation the 
conductivity vanishes, but particle creation at reheating 
quickly turns the Universe into a very good conductor 
at the end of inflation. This short circuits the electric 



3 



field but does not affect the large-scale magnetic field. 
For a close to de Sitter exponentially expanding epoch 
of inflation, on a scale of a thousandth of a Hubble ra- 
dius now (i.e., a few Mpc now), the r.m.s. magnetic field 
strength now is of order 10~ 59 G and very insignificant 
[ill ]. In this case the power spectrum of the magnetic 
field is |i?(fc)| 2 oc fc/a 4 , where k is the coordinate wave 
number, and the evolution with scale factor a is as ex- 
pected from flux conservation. This tiny value is a conse- 
quence of the conformal invariance of the abelian vector 
field lagrangian density and the background spacctimc. 
To generate a large enough magnetic field requires using 
a model in which conformal invariance is broken [Tol - [l2| . 

A simple way to break conformal invariance is to cou- 
ple the inflaton scalar field <j> to the abelian vector field 
through a dilaton-like coupling, generalizing the Maxwell 
lagrangian density to e a ^ F^F* 1 " / e 2 (a is a parameter 
and the exponential form of the dilation-like coupling was 
chosen for simplicity) [HI, [HJ . At the end of inflation, 
during reheating, freezes at its vacuum expectation 
value and the usual abelian vector field lagrangian den- 
sity is recovered. During inflation however, the model be- 
haves as though it has a varying abelian vector field cou- 
pling constant. This modifies the scale factor and wave 
number dependence of the generated magnetic field. For 
power-law inflation 29], power is shifted to the infrared 
and the magnetic field energy density redshifts slower 
than the usual uncoupled abelian vector field case. 

Depending on the value of a, this can result in a 
large enough seed magnetic field to explain the observed 
galactic magnetic fields. The case of greatest interest is 
near the limit of de Sitter exponential expansion infla- 
tion, which results in a close to scale-invariant Harrison- 
Peebles- Yu-Zeldovich spectrum of energy density pertur- 
bations (consistent with the observational indications, 
see, e.g., Ref. [30() with a current epoch magnetic field 
scale-invariant power spectrum \B(k)\ 2 cx fc~ 3 /a 4 , or 
(S 2 (r)) 1 / 2 cx r, and r.m.s. amplitude of order a few nG 
on a scale of a thousandth of a Hubble radius [111 EH ■ 
These computations have been carefully checked, con- 
firmed, and extended [3lj ]. This is an observationally 
viable model. 

In this model, there are regions in model parameter 
space where the vector field fluctuations during infla- 
tion are large enough to invalidate the linear perturba- 
tion assumption, but this does not occur in the region 
of parameter space where a strong-enou gh c urrent-epoch 
large-scale magnetic field is generated [ill . [l2l [H, |33j] . 
That is, backreaction is not a significant issue for this 
classically consistent model. While the abelian vector 
field coupling becomes large during inflation, this is not 
important for the phenomenological, effective classical 
model 0, [111, [H, H3] . Instead, much like the case of 
the "standard" ACDM cosmological model, it is of great 
interest to try to find a more fundamental, quantum- 
mechanically consistent, model that can give rise to this 
classical, observationally-successful, inflationary magne- 
togenesis model. 



An extension of this model, with two scalar fields in- 
stead of one, is also viable (35[. Here, in addition to the 
usual scalar inflaton field, a new scalar field, coupled as a 
dilaton to the abelian vector field, is introduced. Alter- 
natively, as discussed next, the second field is taken to be 
a pseudoscalar, coupling like an axion to the abelian vec- 
tor field. For other inflationary magnetogenesis variants 
see Refs. [36j . 

B. Helical magnetic field 

The inverse cascade scenario, a much discussed mag- 
netic field amplification mechanism, requires non-zero 
magnetic helicity, if it is to be able to transfer magnetic 
power from small scales to large scales and so amplify the 
large-scale magnetic field [37j ■ It is therefore of interest 
to consider inflationary magnetogenesis models that can 
also generate magnetic helicity. 

In a Friedmann-Lemaitre- Robertson- Walker cosmolog- 
ical model, magnetic helicity is [38] 

H B {t) = J d 3 xA • V x A, (1) 

where the vector potential A is related to the mag- 
netic field through a 2 B — V x A. Magnetic helicity 
is usually associated with a nontrivial configuration of 
the magnetic field where the magnetic flux tubes are 
twisted and/or linked [26| . Magnetic helicity resembles 
a Chern-Simons term, for its time variation is given by 
H B (t 2 ) - H B (t x ) = if^d^xF^F^ [38], where F^ is 
the vector field strength tensor dual. This similarity is 
strengthened by the fact that a helical magnetic field may 
be identified as the projection of a non-abelian gauge 
field configuration carrying a nonvanishing Chern-Simons 
number [39j, and that the magnetic helicity coincides 
with this winding number. 

Moreover, helical magnetic fields possess a num- 
ber asymmetry between left-handed and right-handed 
abelian vector field helicity states [H)]. To see this it is 
useful to go to Fourier space and introduce the orthonor- 
mal helicity basis {e + ,e_,e3}, with e± = (ei ±ie2)/V2 
and e3 = k//c, such that {ei, B2, £3} form a right-handed 
orthonormal basis. Decomposing the magnetic field in 
polarization states as B(k) = B + e + + £?_e_, the mag- 
netic helicity reads 

Hs{t) = ^ Jdkk{\B + f-\B_\ 2 ). (2) 

It is clear, accordingly, that a helical magnetic field pos- 
sesses a nonzero difference between the number of left- 
and right-handed abelian vector potential polarization 
states. 

After reheating, the Universe is a good conductor and 
so if a helical magnetic field is created during inflation its 
magnetic helicity survives to the present. Since magnetic 
helicity is odd under discrete P and CP transformations, 
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a cosmological helical magnetic field would be a signature 
of macroscopic P and CP violation. 

Models for generating helical magnetic fields in the 
early Universe exist in the literature [IcJ I39M44} . In or- 
der to generate a helical magnetic field during inflation, 
one can add an interaction term of the form I((f>)F^ v F^ v 
to the abelian vector potential lagrangian density. Here 
I(<p) is a pseudoscalar function of some dynamical or 
background field <j). 

A coupling /(</>) oc 0, between the abelian vector 
potential and the axion field (f> was studied in Refs. 
fill liol . l4lj| . while in Ref. the pseudoscalar <fi was 
assumed to drive inflation. In these cases, because of the 
extra derivative compared to the e°"^ F^ V F^ / e 2 term, 
magnetic power is concentrated on smaller scales, with 
insignificant magnetic power on cosmological scales of in- 
terest. 

In Ref. [43[ , to produce larger-scale magnetic field 
power, the function / was taken to be a time-dependent 
function peaked at long wavelengths (this particular cou- 
pling to the abelian vector potential could be realized by 
a tachyonic massive pseudoscalar field or a massless pseu- 
doscalar field non-minimally coupled to gravity). It was 
shown that, depending on the strength of the coupling, a 
maximally helical field with a scale-invariant spectrum, 
|_B(A;)| 2 oc k~ 3 , could be produced as an excitation of the 
vacuum during inflation. Although it can be quite large, 
large enough to act as a seed for the magnetic fields we 
observe today in galaxies and galaxy clusters, as before, 
its backreaction on the development of inflation is com- 
pletely negligible. 



III. COUPLING OF MAGNETIC FIELD AND 
TURBULENT MOTIONS 

Because of conformal invariance, the usual flat space- 
time relativistic MHD equations are identical to the 
MHD equations in an expanding Universe with zero spa- 
tial curvature when physical quantities are replaced by 
their co-moving counterparts and conformal time r\ is 
used in place of physical time 45j. Based on this fact, 



we perform direct numerical simulations of MHD tur- 
bulence in an expanding Universe using the usual flat 
spacetime MHD equations with a relativistic equation of 
state. Note that our simulations are based on the rela- 
tivistic equations even when studying evolution of turbu- 
lence with non-relativistic bulk velocities. 

Let us briefly describe the primordial magnetic field 
coupling to the fluid. As noted above, the typical char- 
acteristic length scale of this coupling is the phase tran- 
sition bubble size, Ao = 7A//, where 7 < 1 is a parameter 
connected with the number of bubbles, N, within a Hub- 



ble radius, i.e. 
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N 3 , and the Hubble radius 



X H = 5.8 x 10" 



Mpc 



100 GeV^ ^100 



1/6 

Kir) ■ < 3 » 



and QCD phase transitions. The phase transitions are 
characterized by different maximal correlation lengths 
due to the difference in the 7 parameter (that is equal to 
0.01 and 0.15 for electroweak and QCD phase transitions, 
respectively) and in the Hubble radius; see Eq. ([3]) (that 
is equal to 0.006 /ipc and 5.5 pc for electroweak and QCD 
phase transitions, respectively). As we noted above, the 
initial value of the Alfven velocity depends weakly on 
(va oc 1 ^ 6 ), which also makes a difference between 
the electroweak and QCD phase transitions of the or- 
der of unity (1.37), which we discard. The conductivity 
of the Universe is high enough during the phase transi- 
tions: of course the physical characteristics of the plasma 
(such as viscosity and conductivity) depend on temper- 
ature and vary from 100 GeV (electroweak phase transi- 
tion) to 0.15 eV (QCD phase transition). On the other 
hand, in our simulations we assume that, for the goal of 
our study, the main difference between the electroweak 
and QCD phase transitions consists only in a difference 
of the initial correlation length. Our assumption can be 
justified as follows: (i) the forcing amplitude at the ini- 
tial moment weakly depends on the relativistic degrees 
of freedom, oc g+ 1 , (ii) the physical conditions of the 
Universe are different, but new simulations for a wide 
range of the Prandtl numbers show that the growth of 
the correlation length does not depend on the value of 
the Reynolds number, as long as it is large enough |46| . 

We consider two types of forcings, irrotational and vor- 
tical, and for both we assume that the forcing scale coin- 
cides with the bubble size, Ao. The corresponding wave 
number is fco = 2tt/\q. Energy is being dissipated both 
viscously (characterized by the viscosity v) and ohmically 
(characterized by the magnetic diffusivity A, which is in- 
versely proportional to the conductivity). Throughout 
this work we assume that the magnetic Prandtl number 
is vj\ = 1. The smaller the value of v, the more ex- 
tended the turbulent cascade, and the larger the minimal 
mesh resolution required. For all runs presented here we 
use 512 3 mesh points. The largest value of the Reynolds 
number based on the wavenumber fco, Re = u rms /vko, is 
around 200 in all cases (here u rms is the r.m.s. velocity). 

There are many MHD simulations that use vortical 
forcing [23[ but fewer that use irrotational forcing [47| . 
However, we are unaware of any simulation with an initial 
scale-invariant magnetic field spectrum, |_B(fc)| 2 oc k~ 3 , 
of appreciable strength. 3 We generate the correspond- 
ing magnetic vector potential for such a field in Fourier 
space using modes with random phases and suitable am- 
plitudes, as was done in Ref. [48[ for the initial ve- 
locity field using the Pencil Code [49], which is also 



As we noted above we will consider both electroweak 



3 We require that the magnetic field satisfies the upper limits 
from the current CMB and large scale structure observations, [T^l 
(that automatically satisfy the BBN limits), so to be order of few 
nG. Such a magnetic field can be generated during inflation ll|. 
For helical magnetic fields we assume that they are of maximal 
helicity 
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used here. We present magnetic and kinetic energy 
spectra normalized such that j Ey±{k) dk = (B 2 )/2 and 
J E K (k)dk — (u 2 )/2, respectively. 

In all cases we keep the forcing amplitude at a fixed 
level such that the resulting r.m.s. velocity is around 
0.1c s , where the constant c s is the isothermal sound 
speed. The rms value of initial magnetic field strength, 
-Brms, is chosen such that Brms/c s is between 0.05 and 
0.09. 

In the following we describe the types of forcing applied 
to the system. 



A. Irrotational (potential) forcing 
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We consider two types of irrotational forcings. In both 
cases the forcing function is written as f(x,t) — V</>, 
where 0(a;,i) is a random scalar function. In the first 
case we model irrotational forcing in the form of spherical 
expansion waves that reflect the dynamics of phase tran- 
sition bubbles. The forcing and some results for cases 
with no or weak magnetic fields are described in Refs. 
47, 50]. The nondimensional radius of the expansion 
waves is k\R = 0.133, corresponding to a nominal forc- 
ing wave number fco/fci = 15. Here, k\ = 2n/L is the 
minimal wave number in our computational cube of size 
L 3 . In the second case we use random plane waves with 
a wave number k$/k\ = 30. 



B. Vortical forcing 

Vortical forcing is accomplished by generating a ran- 
dom vector potential ip such that / = V x tp. Normally, 
xp is non-helical, but in some cases we have arranged 
such that it consists of positively polarized waves 
Like in the second case with potential forcing, the forc- 
ing wavenumber is here chosen to be fco/fci = 30. 



IV. RESULTS 



A. Blob-like forcing 



Blob-like irrotational [47J and plane wave [5l| forcings 
have been used in the past to investigate the production 
of vorticity through the viscous force, or through the in- 
teraction of irrotationally- forced flow with global rotation 
or shear with an isothermal equation of state, or through 
the baroclinic term with the more general perfect gas 
equation of state [50j . In the present case, again using 
an isothermal gas, the Lorentz force associated with the 
initial magnetic field also produces vorticity [52j . 

In all cases investigated here, the magnetic energy den- 
sity decays, so there is either no dynamo action, or the 
initial magnetic field is still too strong for dynamo action 
to occur because of excessive backreaction on the flow via 



FIG. 1: Spectra of magnetic and kinetic energy as well as 
the kinetic energy of the vortical component at the end of the 
run for vortical blob-like forcing. Re « 250 and Brma/c B 0.05. 
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FIG. 2: Dependence of normalized r.m.s. vorticity on nor- 
malized magnetic energy density during the decay for the run 
shown in Fig. [T] for Re « 250. 



the Lorentz force. As the magnetic field decays, the level 
of vorticity also decreases. In fact, at the end of the sim- 
ulation, the spectrum of the vortical part of the kinetic 
energy follows closely that of the magnetic energy; see 
Fig- HI The r.m.s. vorticity turns out to be approximately 
proportional to the magnetic energy density. Quantita- 
tively, we find w rms /u rms fc « 0.1(£ m /£k) ' 85 ; see Fig. [2] 
Here w rms is the r.m.s. vorticity and £m and £k are the 
magnetic and kinetic energy densities. (We note that the 
exponent 0.85 is probably not robust; in another simu- 
lation we have found a somewhat larger exponent w 1.2, 
for example.) 

During the rising phase of the kinetic energy, the ki- 
netic energy spectrum is approximately proportional to 
k 2 . At small scales, we have an approximate k~ 2 spec- 
trum, in agreement with earlier studies [13, Here 
the spectral magnetic energy exceeds the total spectral 
kinetic energy, and the vortical part of the kinetic energy 
spectrum is only slightly below that of the total kinetic 
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FIG. 3: Magnetic (solid lines) and kinetic (dashed lines) 
energy spectra in regular time intervals. Re « 170. The 
magnetic and kinetic spectra at the last time are addition- 
ally marked in red and blue, respectively. Re ~ 180 and 
B&/C.0.09. 



energy spectrum. 



B. Monochromatic vortical forcing 

Next we consider the case of monochromatic vortical 
forcing. In the simulations presented here, the spectrum 
of kinetic energy is below that of the magnetic energy, but 
it still shows a k 2 behavior at small wave numbers and in- 
termediate times. At earlier times the spectrum is closer 
to a linearly increasing one; see Fig. [3] At intermediate 
times there is a characteristic decline of magnetic energy 
at intermediate wave numbers (k/ko ~ 0.3), which then 
leads to a similar decline of the kinetic energy at these 
wave numbers. 

Visualizations of the B x component of the magnetic 
field and of the logarithmic density, In p, are given in 
Fig. H] for early and late times. The initial fc _1 magnetic 
energy spectrum manifests itself in the form of large ran- 
dom patches which also give an imprint on In p. However, 
as the forcing proceeds, small-scale structures of scale Ao 
become visible in lnp as well as in B x . This is quite 
different in the case of irrotational forcing, of which we 
describe the results from a plane-wave forcing formula- 
tion next. 



C. Irrotational plane- wave forcing 

The irrotational forcing considered in Sec. lIV Al is more 
realistic for applications to phase transition bubbles than 
is plane wave forcing, but to ease the comparison with 
the vortical forcing considered in Sec. IIVBI we now con- 
sider the case of plane-wave irrotational forcing. In this 
case the spectra show initially the same behavior as in 
the vortical case. Even at late times there are similari- 




FIG. 4: Visualization of the B x component and In p for 
the run shown in Fig. [3] Note that at late times the small- 
scale velocity variations cause a corresponding imprint on the 
magnetic field structure. iirms/c s = 0.05 and Re « 180. 
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FIG. 5: Similar to Fig.0 but for potential forcing with plane 
waves. Re « 220 and B r ( m S / c s0.09. 



ties, except that the magnetic and kinetic energy spectra 
now lack the characteristic decline that was visible in the 
vortical case at k/ko ~ 0-3; see Fig. [5] This is because the 
irrotational part of the flow does not interact directly and 
sufficiently strongly with the magnetic field. This is also 
clear from visualizations shown in Fig. |6] which demon- 
strate that the small-scale structures of scale Ao leave no 
imprint on the magnetic field. 

In Fig. [7] we compare spectra of magnetic and kinetic 
energy at early and late times for vortical and irrotational 
plane wave forcing. At the early time the spectra are in- 
dependent of the nature of the forcing, so we show here, 
in the upper panel, only the vortical case at t = 0.5/c s fci. 
At a later time (t = 10/c s ki), the kinetic and magnetic 




FIG. 6: Visualization of the B x component and In p for the 
run shown in Fig. [5] Potential forcing (plane waves). Note 
that at late times the small-scale velocity field has hardly any 
impact on the magnetic field. Urms/cs = 0.07 and Re « 220. 



energy spectra show a decline at intermediate wave num- 
bers for the vortical case (see middle panel), as already 
discussed above. The result for the helical case is virtu- 
ally indistinguishable. This is because magnetic energy 
tends to drive the transfer of magnetic energy to larger 
scales, but the spectral magnetic energy at those scales is 
now already rather strong. In the irrotational case (see 
bottom panel) the decline at k/k 0.3 does not exist. 
This is because here vorticity is small and only this vor- 
tical part of the flow interacts with the magnetic field 
and leads then to turbulent diffusion of the field through 
mixing. 



D. Comparison with helical forcing 

As we have seen in Fig. [71 there is virtually no dif- 
ference between vortical forcings with and without he- 
licity. This is probably related to the fact that there 
was not enough time for the helical forcing to affect 
magnetic and velocity fields at scales larger than the 
injection scale. This can be seen from Fig. [SJ where 
we compare cases with helical and non-helical forcings 
for magnetic and kinetic helicity spectra, JJm(A;) and 
Hn(k), respectively. These spectra are normalized such 
that jH M (k)dk = (A • B) and / H K (k) dk = (to ■ u), 
where ui — V x u is the vorticity. These spectra are 
shown together with their respective magnetic and ki- 
netic energy spectra, Eyi{k) and Ej<z(k), respectively. 
They satisfy the realizability conditions, Em > ||A;-Hm| 
and Ek > \^H^/k\, respectively [53J. (Note that the def- 
initions of J if m and J Hk relative to those of Em and 
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FIG. 7: Comparison of magnetic (solid lines) and kinetic 
(dashed lines) energy spectra at early and late times for vor- 
tical and irrotational plane-wave forcing. The vortical case 
with helicity is shown in light yellow, and nearly coincides 
with the non- helical case (dashed). At early times the spec- 
tra are independent of the nature of the forcing. 

Ek differ from each other by a k 2 factor, which explains 
the slightly different from of their respective realizability 
conditions.) 

Figure [5] shows that in the case with helical forcing, 
where the kinetic helicity is positive, ^kHM(k) is also 
positive and close to E^ik) for k ~ ko, but about an 
order of magnitude below it for smaller values of k. This 
supports the suggestion that helicity effects are too small 
to modify the energy spectra significantly. 

V. CONCLUSION 

In this paper we have studied the evolution of an 
inflation-generated magnetic field [111 coupled to the 
fluid during cosmological phase transitions. Our formal- 
ism is very general and applies to the electroweak and 
QCD phase transitions. The difference between these 
(and other) phase transitions is encoded in the difference 
in parameters such as the temperature and the number of 
rclativistic degrees of freedom, parameters which deter- 
mine the characteristic length scale of the system under 
consideration (Ao). We consider different types of forcing 
and show that at late times the kinetic energy spectrum 
depends sensitively on the forcing used. 

Our forcing scale is determined by the phase transition 
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FIG. 8: Comparison of normalized magnetic (upper panel) 
and kinetic (lower panel) helicity spectra for cases with he- 
lical (dashed) and non-helical (dotted) forcings. In the case 
with helical forcing (dashed lines), red/blue segments indicate 
positive/negative values of Hm and Hk- For comparison, we 
also show magnetic and kinetic energy spectra in the upper 
and lower panels, respectively. 



bubble size. Within a few turnover times the kinetic en- 
ergy spectrum starts to rise on large scales, generating 
large-scale turbulent motions in the fluid. Even a rapid 
phase transition generates turbulence, which will slowly 
decay on large scales. Phase transition-generated MHD 
turbulence might be relevant for cosmological magneto- 
genesis [|[ . Phase transition turbulence can also generate 
a gravitational wave signal that is potentially detectable 

In contrast to previous studies, the inflation-generated 
magnetic field is not frozen into the cosmic plasma. The 
forcing that we considered here is limited by the dura- 
tion of the phase transition. After the forcing source 
stops to act, both magnetic and kinetic energies start to 
decay freely. The configuration of the magnetic field at 
large scales (outside the phase transition Hubble radius) 
is almost unchanged. At intermediate scales correspond- 
ing to the phase transition bubble size there is a slight 



suppression due to energy conversion into kinetic energy. 
The induced turbulent motions are causal so the spectral 
shape at large scales is given by a white noise spectrum 
Ex(k) oc k 2 [55] ]; the vorticity energy density spectrum 
will be steeper (fc 4 ) due to the additional requirement 
of causality (56|. The presence of magnetic helicity does 
not significantly change the forcing stage. On the other 
hand, the scaling laws in the decay stage are strongly 
affected by the presence of magnetic helicity. The dura- 
tion of the decay stage is much longer than the forcing 
stage. During this stage the correlation length of the ve- 
locity increases with a corresponding decay of the total 
energy density. The magnetic field on super-Hubble ra- 
dius scales is decoupled from the fluid which, in turn, 
stays almost unaffected. 

The main results of our study are: (i) inflation- 
generated magnetic fields are not significantly modified 
on large scales by their coupling to the plasma during 
a cosmological phase transition; (ii) the coupling of the 
magnetic field with the phase transition fluid leads to de- 
viations of the magnetic field spectrum from the initial 
scale-invariant shape on intermediate scales; and, (iii) 
there is the possibility of having large-scale correlated 
turbulent motions in the early Universe which, eventu- 
ally, could affect the development of large-scale structure 
formation at late times, and in particular cluster physics 
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